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A theory of nonlinear dynamics of mixed Abrikosov vortices with Josephson cores (AJ vortices) 
on low-angle grain boundaries (GB) in superconductors is proposed. As the misorientation angle 
i? increases, vortices on low-angle GBs evolve from the Abrikosov vortices with normal cores to 
intermediate AJ vortices with Josephson cores, whose length I along GB is smaller that the London 
penetration depth A, but larger than the coherence length f. Dynamics and pinning of the AJ 
vortex structures determine the in-field current transport through GB and the microwave response 
of polycrystal in the crucial misorientation range i? < 20 — 30° of the exponential drop of the local 
critical current density Jb($) through GB. An exact solution for an overdamped periodic AJ vortex 
structure driven along GB by an arbitrary time dependent transport current in a dc magnetic field 
H > Hd is obtained. It is shown that the dynamics of the AJ vortex chain is parameterized by 
solutions of two coupled first order nonlinear differential equations which describe self-consistently 
the time dependence of the vortex velocity and the AJ core length. Exact formulas for the dc flux flow 
resistivity Rf(H), and the nonlinear voltage-current characteristics are obtained. Dynamics of the 
AJ vortex chain driven by superimposed ac and dc currents is considered, and general expressions for 
a linear complex resistivity R(u>) and dissipation of the ac field are obtained. A flux flow resonance 
is shown to occur at large dc vortex velocities v for which the imaginary part of R(uj) has peaks at 
the "washboard" ac frequency u>o = 2nv/a, where a is the inter vortex spacing. This resonance can 
cause peaks and portions with negative differential conductivity on the averaged dc voltage-current 
(V-I) characteristics. Ac currents of large amplitude cause generation of higher voltage harmonics 
and phase locking effects which manifest themselves in steps on the averaged dc I-V curves at the 
Josephson voltages, nhui/2e with n — 1.2, .... 

PACS numbers: PACES numbers: 74.20.De, 74. 20. Hi, 74.60.-w 



I. INTRODUCTION 



Mechanisms of current transport through grain bound- 
aries (GB) in high-temperature superconductors (HTS) 
have attracted much attention because, they reveal the 
d-wave symmetry of the HTS pairinga and determine 
the current-carrying capability of HTS materialsa. Un- 
like low-T c superconductors, GBs in HTS exhibit weak 
link behavior due to the exponential drop of the local 
critical current density of a GB, Jb = Jo exp(— $/-do), 
as the misorientation angle i? between the neighboring 
crystallites increases above $o — 5 — 6°. The strong de- 
pendence of Jf,($) on i? makes high-angle GBs crucial 
current-limiting defects in HTS polycrYstalsa. Since pi- 
oneering experiments of the IBM groupu, much progress 
has been made in understanding the multiscale mi- 
crostructtw;c,|-pf GBs and its effect on their weak-link 
behaviomfflQ, but primarily in the absence of a strong 
magnetic field H. Detailed atomic structure of GBs re- 
vealed by high-resolution electron microscopy have been 
used to determine local under doped states of GB, defect- 
induced suppression of superconducting properties at the 
nano&cale-and controlled increase of Jf, by overdoping of 
GBuffiiM. Recent models have also pointed out the im- 
portance of charging and strain effects which drive the 
HTS state at rGJBi toward the metal insulator transition 
as $ increasesll3Ej. 



At the same time, little is known about vortices on 
low-angle GBs, although it is the dynamics and pin- 
ning of the GB vortices, which mostly limit optical cur- 
rents of HTS polycrystal in a magnetic field 
termine a microwave response of HTS filmsEf 
Generally, pinning of vortices along GBs is weaker 
than in the grains, so GBs form a natural percolat- 
ing network for nrftfenaniial .mntimi-af vortices through 
a superconductoSemBHe'BEf. Such percolat- 
ing networks are not only characteristic of polycrys- 
tals, but represent a rather generic feature of vortex dy- 
namics and pinning in superconductors. For instance, 
networks of easy flow vortex channels have been sug=. 
gested by Kramer in early shear models of flux pinningc3 
and later discovered in molecular dynamics simulation 
of vortices in random pinning potentiaEj'EH, numerical 
simulations of time-dependent Ginzburg-Landau equa- 
tions that [describe moving vortex structure near twin 
boundaries^! and observed by decoratioi£3 and Lorentz 
microscopyBJ. Matching effects in dynamics and pinning 
of mesoscopic vortex flow channels in artificial thin film 
superconducting structures ipse been extensively stud- 
ied by Kes and co-workersE-H. Many observable fea- 
tures of global current-voltage characteristics, magneti- 
zation, rf response and flux creep of HTS polycrystal 
may be due to dynamics and pinning of vortices in easy- 
flow channels, rather than stronger pinned vortices in 
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the grainsii0S'@ii. By contrast, GBs in low-T c ma- 
terials do not bloc^mncroscopic currents, but can en- 
hance flux pinningcic3 and play the role of "hidden" 
weak links which strongly affect the vortex mass and 
viscosityll3 crucial for, transport and microwave response 
of superconductorsEJIl3. 

The behavior of vortices in easy flow channels on GBs 
in polycrystal is mostly determined by the structure of 
vortex cores which depends on the local depairing current 
density Jb through a GB at the nanoscale of few current 
channels between the dislocations. The extreme sensi- 
tivity of J(,(i?) to the misorientation angle # makes GBs 
in HTS a unique tool to trace a fundamental transition 
between the Abrikosov (A) and Josephson (J) vortices. 
As i? increases, J&(i?) rapidly decreases, from the bulk de- 
pairing current density Jd at i? <C i?o down to much lower 
values Jb -C Jd at d ^> i)q. In turn, vortices on a GB 
evolve from the A vortices with normal cores at i? -C i?o 
to intermediate Abrikosov vortices with Josephson cores 
(AJ vortices)t£l and then to the Josephson (J) vortices at 
higher There is no order parameter suppression in the 
AJ core, which is a phase kink whose length I along GB 
is greater than the coherence length £, but shorter than 
the scale of circulating screening currents set by the Lon- 
don penetration depth A. As i? further increases, the A J 
vortices turn into J vortices in which both the Joseph- 
son currents and the magnetic field H(x) vary on the 
same scale, along GB set by the Josephson penetration 
depth AjEJ'E-i This continuous A to AJ vortex transi- 
tion occurs as the spacing between GB dislocation cores 
becomes shorter than £, giving rise to a suppression of 
the amplitude A of the order parameter-^/ = Aexp(iip) 
in current channels between dislocationsEj. Thus, a low- 
angle GB behaves as a high- Jb superconducting-normal- 
superconducting (SNS) Josephson contact, for which the 
Josephson cores of the AJ vortices do not cause pair 
breaking effects responsible for the suppression of A in 
the normal A cores. Such contacts are described by inte- 
gral e ouat jpji.s | p f || a- 1 pQnLQcal Josephson electrodynamics 
(NJE)m^'BBElBBEl which account for the varia- 
tions of phase difference 0(x) = fi ~ (f2 along a GB on 
any length scale greater than £. If 8(x) varies slowly on 
the scales ~ A, the NJE equations reduce to the ps.ua.1 
sine-Gordon equation for long Josephson junctionsc3cj. 
The key difference of the nonlocal approach from the lo- 
cal sine-Gordon theory is that the NJE equations can 
describe the AJ vortex core in the region of parame- 
ters where Jb > Jd£/A, and 9(x) varies on the scale 
I = Aj/A ~ £,Jd/Jb much shorter than the decay length 
A of the circulating supercurrents. 

The importance of the Josephspa-nonlocality for thin 
films has been recognized long agoOEJ. Because the pen- 
etration depth A = 2X 2 /d increases as the film thickness 
d decreases, the nonlocality condition Aj < A can be 
fulfilled even for comparatively low- Jt, junctions. Inde- 
pendently, the NJE approach was developed for bulk su- 
perconductors to describe mixed A J vortices on high- Jb 
" hidden weak links" , such as low-angle GBs in HTS and 



thin a-Ti ribbons in NbTiEl A nonlocal generalization of 
the sine- Gordon equation was also considered in Ref. [f(| 
It turns out that, in the strong nonlocality limit, the NJE 
equation reduces to the well-studied Peierls equation p£ 
dislocation theory, thus, the A J single vortex solutiontp 
is similar to that for the core of an edge diskie&tioti p. 
The NJE equations have other exact solutional2rc3HE3'E3 
for static and dynamic AJ vortex structures. Recently, 
the existence of AJ vortices in low-angle YBa 2 CuzOj 
bicrystals was proven by transport measurements, using 
an exaqt-i expression for the flux flow resistivity of AJ 
vorticesEj. The good agreement between the theory and 
experiment made it possible to extract the core length 
1{T) and the intrinsic depairing current density Jf,(i?) of a 
GB on a nanoscale of few dislocation spacings. The tem- 
perature dependence of Jb oc (T c — T) 2 extracted from 
these measurements does indicate the SNS coupling on 
GBs in HTS, in agreement with the model of Rcf. [fj. 

This paper presents a theory of a nonlinear flux flow 
of AJ vortices driven by dc and ac currents in a mag- 
netic field. An exactly solvable model, that describes 
a chain of AJ vortices moving along a GB through the 
strongly pinned A vortex lattice in the grains, is pro- 
posed. This model of the overdamped vortex dynam- 
ics describes sclf-consistcntly both nonlinear dissipative 
processes in the AJ vortex cores and magnetic interac- 
tion between AJ vortices, showing how the distributions 
of circulating superconducting and quasiparticle currents 
change as a function of the vortex velocity. The AJ 
vortices exhibit many characteristic features of the dy- 
namics of thej-neriodic A vortex lattice, for example, vis- 
cous flux flowE.3, static and dynami c matc hing effects and 
Josephson-like voltage oscillationsE3tiMa. At the same 
time, the AJ vortices can also exhibit effects character- 
istic of the dynamics of short Josephson contacts in ac 
field, or J vortices in long-Josephson junctions, for esapir 
ple, flux flow resonanceEZI, phase locking in ac fielcOEa, 
etc. Pronounced resonance effects occur if the strongly 
overdamped AJ structure is driven by superimposed ac 
and dc currents J(t) = J + J a cos uit at Jo > Jb- It is 
shown that all this rich AJ vortex dynamics is described 
by two coupled first order nonlinear ordinary differential 
equations for the vortex velocity and the core size, for 
any time-dependent J(t). These equations have the form 
of a complex resistively-shunted junction (RS J) equation 
for a short Josephson contact. 

The paper is organized as follows. In Sec. II a qual- 
itative description of length scales of vortices on a GB 
and their evolution with d is given. In Sec. Ill the NJE 
equations and an exact solution that describes a chain 
of AJ vortices driven by an arbitrary ac current in a dc 
magnetic field are presented. In Sec. IV, the nonlinear 
V — J characteristics and the flux flow resistivity of a 
GB in a magnetic field are calculated. In Sec. V a linear 
complex resistivity and rf dissipation are calculated for 
a chain of A J vortices driven by superimposed ac and dc 
currents. A flux flow resonance is predicted. Sec. VI 
is devoted to nonlinear effects caused by superimposed 
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ac and dc currents, in particular the averaged dc V — J 
characteristics in the presence of an ac signal, generation 
of higher harmonics and phase locking effects. Sec. VII 
concludes with a discussion of the obtained results. 



II. VORTEX LENGTH SCALES ON GRAIN 
BOUNDARIES IN A MAGNETIC FIELD. 

The results of this paper are. independent of the de- 
tailed atomic structure of GBS'QEl, so we consider a sim- 
plest planar [001] tilt GB between two crystallites mis- 
oriented by the angle $. Such low-angle GB can be re- 
garded as a periodic chain of edge dislocations spaced 
by do = 6/2 sin(i?/2), where b is the Burgers vectorO. 
Because of the proximity of the HTS state to the anti- 
ferromagnetic metal-insulator transition, regions of size 
~ b near dislocation cores are driven into insulating state 
by local nonstoichiometry, strains and charging effectstii 
As $ increases, the spacing do(^) decreases, becoming 
smaller than the zero-T coherence length £o and the in- 
plane Debye screening length at the angle $o ~ 4 — 6°. 
For t3 > t^o, the proximity effect, strain and charge cou- 
pling cause suppression of A between dislocation cores, 
which becomes more pronounced as increases. In this 
model a GB thus behaves as a SNS Josephson contact, 
whose critical current density </{>($) decreases nearly ex- 
ponentially as increases with </{,($) ~ J a at ■& < i9 - 

A magnetic field H above the lower critical field H c i 
produces A vortices in the grain, and vortices of different 
character on the GB, depending on the ratio £/d . For 
9 <C 9o, the GB vortices are-A vortices with normal cores 
pinned by GB dislocationsEHl. As i9 further increases, a 
GB exhibits a continuous transition from metallic to tun- 
neling behavior above $ > $o, similar to high-Jf, SNS 
Josephson junctional for which the normal core of A vor- 
tices disappears if A on the junction drops below a criti- 
cal valuea. For vortices on a GB, the normal core disap- 
pears quite naturally because of the exponential decrease 
of Jb(^). Since vortex currents must cross the GB which 
can only sustain Jb much smaller than the depairing cur- 
rent density Jd, the modulus A of the order parameter 
is unaffected by vortex cores. As a result, the normal A 
core turns into a Josephson core in which the phase dif- 
ference 9(x, t) on GB varies by 2ir over the length I along 
GB, but the amplitude A is independent of x. The phase 
core length I ~ £,Jd/ Jb is greater then £, but smaller then 
the London penetration depth A, if Jb > Jd/^, jjzhere 
k = A/£ ~ 10 2 is the Ginzburg-Landau parameters. As 
$ increases, the core length l(-d) ~ £,Jd/Jb{i}) increases, 
so the GB vortices evolve from A vortices for i9 <C i?o 
to mixed A J vortices at Jd < Jb{i9) < Jd/n. The A J 
vortices turn into J vortices at higher angles, for which I 
becomes greater than A if 

i?>i? J ~i?oln(A/0, (1) 

where A is the magnetic screening length. In bulk sam- 
ples A is the London penetration depth A, while a thin 



film of thickness d < A the magnetic screening length is 
A = 2X 2 /d. For bulk samples, Eq. (g) y ields 6 J - 230 
if #0 = 5°, k = 100, and Jo = Jd- For a film with 
d <C A, the AJ region expands considerably, for exam- 
ple, if d = 0.1A, then <&j ~ 38°. Therefore, the AJ 
vortices exist in a rather wide range of misorientations, 
$0 < 1? < $j ~ 22 — 40° which comprises the crucial 
region of the exponential drop of Jf,($). In this region 
the in-field current transport through GBs in HTS is de- 
termined by dynamics and pinning of A J vortices. 

Due to the lack of the order parameter suppression 
in the AJ cores, AJ vortices can be described by the 
NJE theory which regards the GB as a Josephson con- 
tact whose Jt(i9) can be tuned in a very broad range by 
varying the misorientation angle Once the AJ core size 
I = £,Jd/Jb exceeds the scales set by the coherence length, 
the Debye screening length and the dislocation spacing, 
the structure of AJ vortices is entirely determined by 
the electrodynamics of currents circulating around the 
cores, regardless of the atomic structure of GB, pairing 
mechanisms and the symmetry of the order parameter. 
The current distribution in the AJ vortex is described by 
the universal Josephson and London equations, while the 
GB imposes a boundary condition of current continuity 
for the sum of the Josephson, quasiparticle and displace- 
ments current densities J = Jb sin 6{x, t) + V/R + CV 
crossing the GB, where V is the voltage on a GB, R and 
C are the quasiparticle resistance and capacitance per 
unit area of a GB, respectively. 

The structure of a vortex on a GB is determined by the 
ratio Kb — X/l reminiscent of the GL parameter k = A/£ 
for the A vortices. Here Kb — (A/Aj) 2 ~ KJb/Jd for 
a bulk sample, Kb — dX 3 /X 2 d for a film of thickness 
d < X. The case Kb <C 1 corresponds to the local relation 
H(x) = <fiod x 8 1 '47rA characteristic of low- J;, high-angle 
GB, for which both 9(x) and H(x) vary on the same 
spatial scale. Such GBs exhibit the J vortkes of length 
Xj described by the sine-Gordon equatioro. However, 
the low-angle GB (especially in thin films) correspond to 
Kb > 1 in which case the relation between 9(x) and H(x) 
is nonlocal, and 9{x) and H(x) vary on essentially dif- 
ferent spatial scales I = £Jd/Jb and A, respectively. This 
gives rise to th|8-|mi xe | d AJ . vj Q¥t4ees-|described by integral 
NJE cquationMyMbliaOSBEl which reduce to the 
sine-Gordon equation if 9{x) varies weakly on the scales 
~ A. The NJE equations provide a universal descrip- 
tion of vortex structures on GB for which all microscopic 
details are hidden in the intrinsic parameters Jb and R 
of a GB. These parameters are very difficult to calcu- 
late, given the present state of the microscopic theory of 
HTS, but they can be extracted from resistive measure- 
ments on HTS bicrystrals with the help of exact NJE 
solutions that describe the flux flow resistance of moving 
AJ vorticesEJ. 

In a strong magnetic field H S> H c \, the A vortex 
spacing a = (</>q / H) 1 / 2 is shorter than A, thus the relation 
between 9(x,t) and H(x,y) is always nonlocal, regardless 
of the value of Jb ■ The A J vortex chain is then has two 
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length scales: the core size I > £ and the inter vortex 
spacing a(H). Both lengths a and / become comparable 
at a characteristic field H ~ <j>o/l 2 ~ (•/&/ Jd) 2 H c2 much 
smaller then the upper critical field H C 2- Thus, unlike 
the A lattices, there is a wide field region Hq < H < 
H c 2 in which the AJ vortex cores overlap, but the bulk 
superconductivity persists. 

The larger core of AJ vortices leads to their weaker 
pinning along a GB, which thus becomes a channel for 
motiop_of AJ vortices between pinned A vortices in the 
grainalj (Fig. 1). This gives rise to an extended lin- 
ear region in the V — I characteristic of a polycrystalline 
HTS-tjaa J; i s .de|tMpated by motion of AJ vortices along 
GBs£3E3cafi3'ESc3c3. Pinning of AJ vortices results from 
interaction of the AJ phase core with structural inhomo- 
geneities of GB and the magnetic interaction of AJ vor- 
tices with more strongly pinned A vortices in the grains. 
If the periods of the AJ vortices and bulk A vortices are 
slightly different, the A J vortex chain breaks into com- 
mensurate domains (in which A and AJ periods coincide) 
separated by domain walls. This behavior, characteristic 
of commensurate-incommensurate transitionsEErEil, was 
observed in molecular-dynamics simulations of A vortices 
in artificial flux flow channelsES for which the width of the 
domain walls (dislocations) considerably exceeds the in- 
ter vortex spacing. In this case the pinning of domain 
walls by the intrinsic Peierls potential is exponentially 
weak, so the dcpinning critical current I g b(H) is most 
likely due to macroscopic variations of superconducting 
properties along GB, for example, facet structures which 
cause significant peaks in I g b(H) if the A_vortex spac- 
ing is commensurate with the facet periodo, strains and 
local nonstoichiometry along GB0, etc. 

At low field, only a single A J vortex row moves along 
GB, while the A vortices in the grains remain pinnedEil. 
At higher field the moving AJ vortices start dragging 
neighboring rows of A vortices in a flux flow channel 
along GB. The field Hi below which only a single A J 
vortex row moves along the GB, can be estimated from 
the condition that the pinning force f m of AJ vortices 
due to their magnetic interaction with fixed A vortices 
equals the intergrain pinning force </>qJ c /c. The pin- 
ning force of AJ vortices is the maximum gradient of the 
magnetic energy f[x) — —<f>od x H(x)/ATT, where H(x) = 
B+AH cos(2irx/a) is the local field produced by the fixed 



A vortex lattice along GB, AH 
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amplitude of the oscillating part of the local field H{x) 
due to the periodicity of the A lattice, and-,u ~ a is the 
spacing of the first A vortex row from GBtEl. Therefore, 



27i-A 2 J c i 2 Attu, 

exp( ). 

a 



<V<Po 



(2) 



For u = a, A = 20001, and J c = 10 5 - 10 6 A/cm 2 , Eq. 
(0) yields Hi ~ 0.1 — 10T. Note that the essential depen- 
dence of the transition field Hi on the bulk J c indicates 
that the region < H < Hi can be significantly widen 
by irradiationj-ssjhich increases J c while weakly affecting 
GB properties!^. In addition, Hi is very sensitive to the 



position u of the first vortex row which can be strongly 
affected by the GB microstructure, facets, long-range T c 
variations due to strains, etc. 

In this paper we neglect the pinning of vortices on 
a GB, assuming that the driving current J(t) is higher 
than the depinning current of the AJ vortices, but lower 
than J c of the A vortices in the grains. This behav- 
ior has been observed on HTS J^p CTst|ak -ia- p jT^ e re ~ 
gion H < Hi of magnetic fielddi^^BBEHOM Un- 
der the ac current, pinning effects in electromagnetic re- 
sponse weaken even more as the ac fcesuency uj exceeds 
a characteristic depinning frequency! 3 "HI 65 !. Furthermore, 
low-angle GBs can be regarded as overdamped Joseph- 
son contacts for which the displacements currents can be 
neglected. Indeed, the overdamped state corresponds to 
uj <C uj c = (CR)^ 1 ~ Air/SpSe, where Sp and 8e are the 
excess resistivity and dielectric susceptibility on GB. Be- 
cause for low-angle GB the dielectric dislocation cores do 
not overlap, 5p ~ p ni and Se ~ 1, where p n is the normal 
state resistivity at T c . Thus, the condition uj <C uj c always 
holds for uj smaller than the superconducting gap (see 
also Refs. ^,^). We also neglect the time dispersion of 
the the GB resistance R and contribution of bulk quasi- 
particles, adopting the simplest, frequency-independent 
R in the framework of a standard RSJ modela. 



III. GENERAL DYNAMIC EQUATIONS. 



The NJE equations for current-driven vortex. Stru c 



tures on an overdamped Josephson contact arel 

n ^ : — — r / e'(u)K 
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I = c4,/167r 2 A 2 J fe , 



du + B Vl (3) 

WO.0 + P, (4) 

/2ncRJ b . (5) 



Here the overdot and the prime denote differentiation 
with respect to time and the coordinate x along GB, 
fl(x,t) — J(x,t)/Jb, J = (c/47r)V x B v is the current 
density across GB induced by bulk vortices, 



B v (x,y) = 



2ttA 2 ^ 



An 



r - iv 



(6) 



where r n is the position of the n-th A vortex, Ko(x) 
is a modified Bessel function, 0o is the flux quantum, 
and c is the speed of light. The first term in the right- 
hand side of Eq. (S) describes the magnetic field pro- 
duced by all currents circulating near GB, and the second 
term represents the contribution of bulk vortices with- 
out GB. To provide the boundary condition for 3(x, y) 
on a GB, the phase difference 9(x,t) must satisfy Eq. 
(Q), which results from the current continuity condition, 
cH'/4n = J 6 sin6» + H6/2eR. Eqs. describe spa- 

tial variation of 8(x, t) and H (x, y, t) on any scale greater 
than £, irrespective of the microscopic mechanisms of 
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current transport through the GB. The only assumption 
Jb "C Jd behind Eqs. (H)-(|§) ensures the lack of the order 
parameter suppression by currents flowing through GB. 
The geometry of the sample manifests itself in the long- 
range asymptotics of the kernel in Eq. ([|) at |r — r'| > A. 
For instance, Eqs. (||)-(^|) correspond to an infinite GB 
in a parallel field. More complicated expressions for the 
kernel K[(x — u)/X] have.been obtained for thin films in 
a perpendicular fieldEJcZlca, and slabs in paralleled and 
perpendicularly fields. 

In the vicinity of the AJ cores, r <C A, and also in 
the high-field limit, H 3> H c i, Eq. (^) acquires a simple 
universal form independent of the sample geometry. This 
universality results from the fact that for H 3> H c \, the 
derivative 9"(u) in Eq. (|J) rapidly oscillates over the 
inter vortex spacing (tp^/H) 1 / 2 <§; A. In this case the 
main contribution to the integral comes from the region 
| a; — u\ < A, where the Bessel function K (x) in Eq. (^) 
can be replaced by its expansion Kq(x) ~ — ln(x) at 
small x. Thus, the equation for 9 becomes 



TU = — 
7T 



' {u)du 



u ~ x 



sin ( 



(7) 



For other geometries, Kq[(x — u ) /A] in Eq. (|) should be 
replaced with the appropriate kernel K[(x — u)/X] which 
always has a logarithmic singularity at x = u, and a 
geometry-dependent nonsingular part, K reg (x,u): 

K{\x-u\/X) = -]n\{x-u)/X\+K reg [{x-u)/X]. (8) 

This general behavior of K[(x — it) /A] is illustrated in 
Appendix A where K[(x — u)/X] for a thin film (d <C A), 
is considered. Because of the rapid oscillations of 0"(u) 
at H ^s> H c i, the main contribution to the integral in 
Eq. (|J) comes from the narrow region around u = x, 
so neither K reg [(x — it)/ A], nor the screening length A 
contribute to Eq. (g). This feature of Eq. ^ reflects 
the physical fact that the distribution of currents near 
the AJ cores is unaffected by the London screening. In- 
deed, the Green function K[(x — u)/X] is proportional to 
the single- vortex London solution H(r), so the difference 
between K [{x — u)/X] for bulk samples and thin films is 
basically the same, as between H(r) for the A vortex and 
the Pearl vortexfEl, respectively. Both vortices have the 
same distributions of currents J(r) near the core, r < A, 
but very different asymptotics of J(r) for r > A. Thus, 
the universal Eq. (m) describes the distributions of the 
phase difference 9{x) in the AJ vortex cores and circulat- 
ing supercurrents on the scales < A away from the cores 
where the London screening is inessential. 

The driving parameter [3 — /3q + 8(3{x) is a sum of the 
constant transport current [3q due to the gradient of the 
A vortex density in the grains, and an oscillating compo- 
nent 5j3{x) due to the discreteness of the A vortex lattice. 
The term 5(3(x) gives rise to a critical current (3 C through 
GB due to pinning of AJ vortices by A vortices in the 
grainsEj. In this paper we consider a rapidly moving AJ 



structure in the flux flow state, (3 (3 C , for which the 
pinning term 8(5{x) <§C 1 can be neglected, and (3{x) be 
replaced by 0o(t). As shown in Appendix B, the non- 
linear Eq. (|7j) then has the following exact solution that 
describes a stable periodic vortex structure: 



= 7T + 7 + 2tan~ 1 [Mtan£;(x - x )/2), 



(9) 



where j(t), M(t) and the vortex velocity v(t) — xq de- 
pend only on t and obey the following equations 



tq + sinh a cos 7 = 
7-7 + sin 7 cosh a = /3o(i), 
krv = — sin 7 sinh a. 



(10) 

(11) 
(12) 



Here h — (kl) 2 is the dimensionless magnetic field, the 
wave vector k — 27r/a defines the period a of the A J 
structure, and 



sinha = 2M/(M 2 - 1). 



(13) 



Using the complex variables z = j+ia and / = f3o + i\/h, 
Eqs. (|lO|)-(p^) can be written in a more compact form 



krv — Imcosz. 



(14) 



The equation for the complex "phase" z(t) has the same 
form as the usual RSJ equation for the phase difference 
on an overdamped point contact. 

Eqs. (p|)-(^2|) were obtained in Appendix B by the 
Hilbert transform, which was used to obtain static pe- 
riodic solutionSj-of the Peierls equation ([?]) in the dis- 
location tteoryO, and then employed to describe AJ 
structuresEa. It is instructive to derive Eqs. (|9|)-(|l^) us- 
ing a more transparent approach, starting from the basic 
London equation 



H - X 2 V 2 H = <p 6(x)6(y)/2<K, 



(15) 



where 9[x) is determined by a particular vortex structure 
on a GB. Since screening does not affect the A J cores and 
current distribution on the scales < A away from GB, Eq. 
( jilt ) reduces to the Laplace equation V 2 H = supple- 
mented by the boundary condition for the tangential and 
normal components of the current density on GB: 

d y H{x,+0) - d y H(x, -0) = -<b 9'(x)/(27rX 2 ), (16) 

J y {x, ±0) = ~H'{x, ±0) = J b sin9+ - J. (17) 

47r 2ncR 

If screening is inessential, H(x, y) becomes a potential 
field, which in some cases can be found directly using 
the theory of analytic functions. For instance, H(x,y) 
produced by the periodic AJ vortex structure in the upper 
half-plane y > 0, is given by the following ansatz 

H = ^LRelnsinla; - x {t) + i(\y\ + y (t))]^ (18) 

which describes the field produced by a chain of ficti- 
tious A vortices displaced by y — —yo away from GB in 
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the lower half-plane (we omit a constant term in H (x, y) 
inessential for \y\ < A). Eq. ( fl8| ) is an exact solution of 
the Laplace equation V 2 H = 0. It turns out that yo(t) 
can be chosen such that H (x, y) also satisfies the bound- 
ary conditions (|l6|) and jl7| ) provided that 9(x,t) is g iven 
by Eq. (§) while 7, a, and x (*) obey Eqs. ©-(111). To 
show that we first observe that 



d x H 



k(j> 



sin £ 



47rA 2 [cosh k(y — y ) — cos ' 
5 iJ - fe 0o sign(y) sinh a 

47rA 2 [cosh k(y — y ) — cos ' 



(19) 
(20) 



where £ = [x — xq (t)]k. In turn, Eqs. @ and Q yield 

„, fcsinha 



sin 6* 



cosh a — cos £ 
d sin £ + fcio sinh a 
cosh a — cos £ 
( 1 — cosh a cos £) sin 7 — sinh a sin £ cos 7 



7' 



(cosh a — cos £) 



, (21) 
, (22) 
(23) 



Now 2/0 can be chosen such that the denominators of Eqs. 
(0)-(|23|) would coincide at y = 0: 



fclto = a = ]n[(M+l)/(M-l) 



(24) 



Eqs. (|2C|) and (|21| ) with fcyo = a automatically satisfy 
the first boundary condition (fHf). Furthermore, substi- 
tuting Eqs. ([l9|)-(p4]) into the second boundary condi- 
tion ( |l7j ) reduces the latter to the form: C\(t) cosfc£ + 
Czit) sinfc£ + G$(t) = 0. The self-consistency conditions 
Ci{t) = are satisfied only if 0(t), a(t), and Xo(t) do 
obey the dynamic equations Jl0|)-(p^|). 

Fig. 1 shows the current streamlines calculated from 
Eq. ( |l8| ) , which has a-clear interpretation similar to that 
of a single AJ vortexEfl Namely, the current streamlines 
described by Eq. ( |l8| ) in the upper half-plane y > coin- 
cide with those produced by a chain of moving fictitious 
A vortices displaced by y — —a(t)/k away from GB. 
Likewise, the current streamlines in the lower half-plane 
y < coincide with those produced by a chain of ficti- 
tious A vortices displaced by y = a(t)/k away from GB. 
The resulting nonsingular field distribution H(x,y) is an 
exact solution for moving A J vortices, where the trans- 
verse displacement j/o determines the AJ core size along 
GB. The time-dependent Cartesian coordinates Xo(i) and 
yo(t) = a(t) /k of these fictitious A chains obey Eqs. (|To|)- 
(O). The transition from AJ to A vortices occurs as 
increases, reaching Jb ~ J^, while yo decreases down to 
2/o ~ £. Likewise, the transition from AJ to J vortices 
occurs as Jh decreases below Jd/n, in which case yo be- 
comes greater than A. 

The set of coupled ordinary differential equations (|l0|)- 
( p^| ) describe the evolution of the AJ phase core length, 
the phase shift 7(i), and the vortex velocity io(t) for any 
time-dependent transport current (3q{£). These equations 
along with Eqs. (Q) and fll8| ) determine distributions of 
the phase difference 0{x,t) and screening currents for an 



interacting moving AJ vortex chain, including nonlinear 
dissipative dynamic states caused by ac and dc driving 
currents of large amplitude, J(t) ~ Jf,. Eqs. (p^|)-(p^) 
considerably simplify in the low- field limit, a — > 00 for 
which Eq. (|^) reduces to a_superposition of independent 
single AJ vortex solutionsEj 



0(x J *)=7(t) 



2 tan 



x (t) 



L(t) 



(25) 



Here k = 2w/a — > 0, a — > 2/M — * 0, so the AJ core 
length L(t) = a/k is independent of the magnetic field, 
while a = Lk^0. If a -> 0, Eqs. ©-© turn int&^he 
following equations that describe a single AJ vortej£3 



T7 + sin7 = P , 
tL + Lcosj = I, 
txq = —L sin 7. 



(26) 
(27) 
(28) 



Unlike Eqs. ([^ -( fLf)), the equation for 7 is decoupled, 
from Eqs. (p7) and (g8|. Thus, in the limit a > L, the 
dynamics of the core length L(t) and the vortex veloc- 
ity Xo is integrable for any given 7(f) which in turn is 
determined by Eq. (|2^) for any time-dependent /3o(t). 
This result is no longer valid if the interaction between 
A J vortices at a finite H is taken into account, when 
both j(t) and a(t) are determined self-consistently by 
Eqs. 

Contribution of each A J vortex to the "staircase" so- 
lution (]25|) gives a 2n phase shift along GB, 4-hus the 
AJ vortex carries exactly one flux quantum (/>otil. Gen- 
erally, the inter vortex spacing a(H) on GB is different 
from the period {4>q/ H) 1 / 2 ot-the bulk A lattice, because 
of the reduced H c i on a GBEl. In the low-field region, 
H ~ H c i, the magnetic induction B is very different from 
H, which causes a significant mismatch in the periods of 
AJ and bulk A vortex lattices. However, for strong fields 
H 3> H c i considered in this paper, the difference between 
B and H is negligible, so the A J spacing a(H) nearly co- 
incides with the bulk one, {(po/H) 1 / 2 as both are fixed 
by the same flux quantization condition. 

The averaged voltage V on a GB produced by the mov- 
ing AJ vortex structure is given by: 



V 



2nca 



6dx. 



(29) 



Substituting Eq. fl22|) into Eq. (|29|), we see that the 
term proportional to a in 9 vanishes after integration, 
while the part of 9 (x — xo) proportional to xo reduces 
to the full derivative —±q9'. Integration of 9' fronuX) 
to a gives 2ir due to the flux quantization conditiorcil, 
9(x + a) - 9(x) = 2ir, thus 



V= ^- kv) - 



(30) 



Here the first term in the parenthesis describes the quasi- 
particle component of V, and the second term results 
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from vortex motion. Eq. (|3C|), along with Eqs. (|10|)- 
(|l2|), determine the nonlinear electromagnetic response 
of the moving AJ vortex chain. These equations can be 
presented in different forms, depending on the way the 
external drive is applied. For example, the solutions of 
Eqs. (|l(])-(|lj) determine the voltage ( |30| ) for any time- 
dependent current J(t). 

Another regime corresponds to the rf response, when it 
is the voltage V(t) (rather than the current density J(t)) 
on a GB which is fixed by an external rf source. In this 
case it is convenient to express the vortex velocity via 
V from Eq. © and subtract Eq. (0) from Eq. (§§). 
Then the equations for a(t), "f(t) and f3(t) in the fixed 
voltage mode take the form 



ra + sinh a cos 7 = Vh, 
T7 + sin 7 sinh a = u, 
(3 — u + e~ a sin 7, 



(31) 
(32) 
(33) 



where u = V/RJb is the dimensionless voltage on a GB. 



IV. NONLINEAR STEADY-STATE FLUX FLOW 

The steady-state propagation of the AJ vortex chain 
driven by the constant current /3q is described by Eqs. 
(0) and (|l0|)-(|l^) with a = 7 = 0, whence tanha = 
—svh//3o, tan 7 = —s, and syh — — sin 7 sinh a. These 
equations give the dimensionless propagation velocity 
s (Po) = v / v o m the form 



give v = v (3o/y/l - ffi 
tex velocity v(J) diverges at J 
core size L 



s 2 = [^(l-Pl + hY + APlh -l-h + $]/2h. (34) 

Here vq = 1/t, h — (2irl/a) 2 is a dimensionless magnetic 
field. The limit h — > corresponds to a single AJ vor- 
tex, for which both Eq. (34) and the steady-state Eqs. 
|) (L cos 7 



/, sin 7 = Pq and tv = — LCOS7), 
The so-obtained single vor- 
Jb, because the A J 
l/y/l — (3 2 expands as /?o increases^]. The 
core expansion aSpthe velocity v(J) increases is charac- 
teristic of both A,k3 and Jo vortices in the overdamped 
limit, unlike the Lorentz contraction of J vortices in the 
underdamped limit. However, as the AJ cores expand, 
they start overlapping, so the interaction between vor- 
tices at J > Jb cannot be neglected even at low fields, 
h -C 1. The importance of the interaction is apparent 
form the exact Eq. (|34"|), which shows that the veloc- 
ity s(/?o) smoothly increases as (3q increases and has no 
singularity for any nonzero h. 

The dc V — J characteristic due to viscous motion of 
AJ vortices follows from Eq. ([50|) in which 7 = 0, and 
v(J) is given by Eq. (j34"[). Hence 



V2 



(1 - I3 2 + hf + Affih -l-h + {3 2 



1/2 



where Vq = RJb- Eq. ( p5| ) can also be written in a more 
transparent form 



J 



V 
R 



1 



H/H Q + (V/V ) 2 



1/2 



(36) 



The V — J curve shown in Fig. 2 has two linear portions: 
a flux flow part, V = RfJ at J<C Jb, and the quasipar- 
ticle Ohmic part, V = RJ for J ^> J&. In the crossover 
region, J ~ Jb, the V — J curve becomes nonlinear, be- 
cause of the A J core expansion as J increases. Since 
the A J vortex core overlap at J > Jb, the GB resistance 
approaches its maximum value R at J 3> Jb- 

For J <C Jb, Eq. @ yields V = RfJ, where R f = 
R\Jhj (1 + h) is the flux flow resistivity of AJ vortices. 
If H > Hoi, then h = {2i:l/a) 2 = H/H , and 



Rf = 



rVh 

VH + H ' 



Hn = 



(37) 



At H <C H Q , Eq. ( pq j describes a chain of A J vor- 
tices whose cores do not overlap. In this case Rf(H) 
is similar to the ID Bardeen-Stephen formulae, Rbs — 
R^J H/ H c2 , except that in Eq. (^) the core structure 
is taken into account exactly (see also Ref. |l9|). For 
H > H ~ (Jb/ Jd) 2 H C 2 ^ H c2 , the AJ cores overlap, 
and Eq. ( |37j ) describes a crossover to a field-independent 
resistance R. This regime has no analogs for A vortices, 
whose normal overlap only at H C 2- 

It is interesting to compare Eq. (pJ) to Rf(H) for 
J vortices at J <C Jb- It is knownEJ that for small 
density of J vortices, Rf(B) = RB/H c u is propor- 
tional to B, similar to the Bardeen-Stephen resistivity 
in which H C 2 is replaced by the Joscphson lower critical 
field H c u = (f>o/TT 2 \\j. A general expression for Rf(H) 
can be written in the following form (see Appendix C): 



Rf(H) = RB(H)/H. 



(38) 



(35) 



Here B(H) is determined by. the equilibrium magnetiza- 
tion the J vortex structureES: 

B = Tr 2 H cl j/4pK(p), H = H clJ E( P )/p, (39) 

where K (p) and E(p) are the complete elliptic integrals 
of the first and the second kinds, respectivelyEJ, and 
< p < 1 is a continuous parameter. As follows from 
Eqs. (|3^) and (]39|), Rf first increases linearly with B at 
B <C H c ij and then approaches the constant value R at 
H ^> Hcij as the J vortices overlap, and B — > H . By 
contrast, the saturation of Rf(H) for the A J vortices oc- 
curs at much higher fields ~ Hq ^> H c \ for which the AJ 
vortex spacing a — ((fio/B) 1 / 2 becomes of order the AJ 
core size I. The dependence of Rf(H) for J vortices looks 
rather different from the simple Eq. ([37]) mostly because 
of the complicated relation (|3^) between B and H at low 
fields H ~ H c ij. In fact, Eq. d37j ) also becomes more 
complicated at low H, of the order of the AJ. lower crit- 
ical field H cU = (</> /47rA 2 )[ln(A/Z) + 0.423]E3, for which 
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H in Eq. ( |37| ) should be replaced-by the corresponding 
B(H) dependence for AJ vorticesc£l. For H » H clb , the 
induction B almost coincides with H, thus Rf(H) for A J 
vortices acquires the universal form ( |37| ) independent of 
demagnetizing effects crucial at H ~ H&. The simplic- 
ity of Eq. ( p7|) is very convenient to extract intrinsic 
properties of GBs at the nanoscale from transport mea- 
surements of Rf(H) in HTS bicrystalsEl 



V. LINEAR AC RESPONSE 

To obtain the dynamic linear resistance R u for a weak 
ac current, J a cos(oji) <C Jb superimposed on the dc cur- 
rent J, we calculate the amplitude of the induced ac 
voltage, V u from Eq. (M) in which ^(t) and v(t) are 
determined by Eqs. (|1C|)-(|12[) with (3 — (3q + (3 a exp(iwi). 
Setting a = a + 5a, 7 = 70 + ^7 and calculating the 
perturbations 8a <C 1 and £7 -C 1 induced by the ac 
current from the linearized Eqs. ( 10)-([l2|) yields the fol- 
lowing general expression for the ac complex resistivity 
Rui = Vuj/Ja (see Appendix D): 



Ruj 
~R 



jj{h + u 2 )/Vh - Q 2 + in(r) + Vh) 



fi 2 
"0 



n 2 + 2irjn 



(40) 



Here u = V/ JbR is the dimensionless dc voltage on a GB, 
£1 = lot is the normalized ac frequency, fio and rj = (3q/s 
are the dimensionless flux flow resonance frequency and 
viscosity, respectively 



■q = y/h + h/(h + u 2 ), 



n 



(41) 



In the fixed current mode, it is convenient to express 
Qq and r\ in terms of (3$: 



V 



tt» = [(l + h-f3 2 ) 2 +Ah(3 2 ] 1 /\ (42) 
I + ft — Pq] 1 ^ 2 / V2 (43) 



[ v /(l + ft-/3 2 ) 2 +4ft/3 2 + l + 



The dependencies of and r\ on Pq are shown in Fig. 
3. For J < Jb, £Iq and 77 practically coincide, but for 
higher J, the frequency Qq becomes much higher than the 
damping constant rj. In the following subsections various 
dynamic regimes described by Eq. (J40|) are considered. 



Small dc vortex velocities 



For v — 0, Eq. (|4(]) reduces to 

Rf-R 



R^, — R 



lUJTf 



1 



(44) 



Here Rf is the dc flux flow resistivity (37), and r/ is the 
flux flow relaxation time constant in the AJ vortex chain: 



T / = 



V 7 ! + H/H Q 



(45) 



The Drude-like frequency dependence of R^, results in 
the following retarded relation between the induced volt- 
age V(t) and the driving current J(t) 



V(t) = RJ(t) + (R f - R) 



* ''-« dt' 

e —j(t') — . (46) 

-00 T / 



For example, after a jump- wise increase of J{t) from to 
Jo, the steady-state flux flow sets in according to 



V[t) = J [R f + (R- R f )e^' T f], t > 



(47) 



and V = for t < 0. The discontinuity in V(t) at t = 
disappears if a time dispersion of R, or a finite capaci- 
tance C of the GB are taken into account. In the latter 
case V(t) sharply increases from to V(t) given by Eq. 
([fT]) during a short time = RC. 

The ac power Q = (\/2)JlRe(Ru), dissipated on a GB 
due to viscous flow of AJ vortices can be obtained from 
Eqs. @), (0), and <M) in the form 



Q = 



RJ 2 a [^h(TTh) + {lot) 2 ] 
2[l + h+(ujT) 2 } 



(48) 



For a fixed frequency, Q(H) monotonically increases with 
H, approaching the quasiparticle limit RJ 2 /2 for H ^> 
Hq, as shown in Fig. 4. In the steady state, lot — > 0, 
the power Q{H) oc \fhj\J\ + h is simply proportional to 
Rf(H). For finite frequencies, Q(H) becomes finite even 
with no vortices (H = 0) due to quasiparticle ac ohmic 
currents through GB. 



B. Flux flow resonance for moving AJ vortices. 

Interaction of the moving AJ chain with the ac field 
can cause a resonance, if lo is close to the real part loq of 
the complex eigenfrequency lo f which corresponds to the 
pole in i? w (f2): 



cof = loq + irjo, 



(49) 



where the flux flow resonance frequency loq and the 
damping constant 770 are given by 



= \J{kv) 2 



no 



H 



1 



1 + (v/v ) 2 



1/2 



(50) 
(51) 



For small dc velocity v, the AJ oscillations are strongly 
overdamped, loq — r/o, so no resonance peaks in i? w (fi) 
occur for J < Jb, as evident from Fig. 3. However 
because the frequency loq increases, while the damping 
constant rjo decrease as v( J) increases, the flux flow res- 
onance emerges at high dc driving currents J > Jb, for 
which loq 770 ■ For a given frequency lo, the resonance 
occurs at the vortex velocity v / , for which 



lj 2 = (kv f ) 2 +rj 2 (v f ). 



(52) 
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If kv 3> t]o, the resonance frequency ujq approaches the 
"washboard " frequencyL3 for which the vortex velocity 
Vf equals the phase velocity of the electromagnetic wave 
oj/k with the wave vector k — 2it /a of the AJ structure. 
Because rio(/9o) has a minimum at J ~ Jf, (see Fig. 3), 
the resonance condition ( j52[) at a given frequency u> can 
be satisfied either for one or two velocities Vf. From Eqs. 
©, and ©, it follows that 



Vf 



Ah 



2h 



1/2 



V . 



(53) 



The solution for Vf exists above the threshold il > 
(4/t) 1 / 4 , where the resonance can occur either at one or 
two different velocities Vf, depending on Q: 



(4/1) 1 / 4 < n < {i + h) 1 ' 2 
n > (i + /i) 1/2 



two Vf (54) 
one Vf. (55) 



At Q = (1 + h) 1 / 2 , the smaller resonance velocities van- 
ishes as Vf oc (VI + h - £1) 1/2 . For > /i 1/4 , the larger 
resonance velocity approaches the material-independent 
"washboard" value Vf = aoj/2ir. 

Near the resonance, u) ~ ujq, at high vortex velocities 



(s 2 > 1, fi « s\/7i, ?y w Vh < ft ), Eq. 



yields 



= 1 



1 



4fi {if] - + ^o) 



(56) 



Eq. ( |56"1 ) describes a resonance line with a Lorentz peak 
in Imi? w = — i?./4s[ft,+ (fi — Sl ) 2 ] , as shown in Fig. 5. For 
s 2 ^> 1, the amplitude of the peak decreases as v(J) and 
H increase. The resonance is most pronounced if J > Jb 
and H <C Hq, while at smaller J or higher H, the peak in 
ImR^^l) disappears as the linewidth r\ becomes of the 
order of the eigenfrequency Qq. 



VI. NONLINEAR RESPONSE 

A. Flux flow resonances on the dc V-J curve 

The flux flow resonance also manifests itself in the aver- 
aged dc voltage V as a function of the dc current density. 
To calculate V(J), it is convenient to use the complex 
representation ( J14| ) for z = r y+ia, taking z = zq+Szq+Sz, 
where zq is the dc solution without the ac field for which 
sin z Q — {3 + i\fh, and the oscillating correction 6z obeys 
the linearized equation 



tSz + cos zqSz = 6(3. 



(57) 



Here 8j3{t) = (3 a cos tot, and 6z is a dc correction to z 
due to the ac field, which is determined by Eq. Jl4| ) 
expanded to quadratic terms in 8z: 



25zqcoszo = (Sz 2 ) sin zq, 



(58) 



where the angular brackets mean time averaging over the 
ac period 2tt/lo. The averaged vortex velocity (v ) is given 



by the second of Eqs. (|lj) expanded to quadratic terms 
in Sz and linear terms in 5z : 



kr(v) = Im ( cos zq 



(62 



2 COS Zq 



(59) 



The second term in the parenthesis describes the correc- 
tion due to the ac field. The value (Sz 2 ) can be calculated 
using the solution of Eq. (|57 



Sz = 



Pa 



, -(cos zo cos ujt + Q sin ujt), (60) 
il z + cos 2 zo 

which yields (Sz 2 ) = B 2 /2[£l 2 +cos 2 zp\. Inserting this ex- 
pression into Eq. (p3) gives the averaged vortex velocity 
in the form 

Vh(s) = -(1 + r u ) sin 70 sinha , (61) 
where the parameter T w quantifies the ac contribution 

j 2 (4 V 2 + n 2 -n 2 ) 



r w = 



(2j b n Q ) 2 [(n 2 



n 2 ) 2 



Ar] 2 n 2 



(62) 



Eq. (|6l| ) f° r ( s ) ls identical to the dc relation, rkv = 
— sin7sinha, if v is replaced by the effective velocity 
v = v/[l + r w (v)]. Therefore, the averaged J — V char- 
acteristics can be presented in the dc form ([36]) if V is 
replaced by the effective voltage V(V) as follows 



1 



1 



1/2 



(63) 



H/H + (V/V ) 2 - 
V = V/[l + r u (V)], (64) 

where T(V) is determined by Eqs. © and ©. The 
J — V characteristics described by Eqs. ( |63|) and (|64| ) 
J a are shown in Fig. 6 for different ac amplitudes J a . 
For stronger ac signal, the J — V curves can exhibit two 
maxima at the resonance voltages Vf — (j)QVf/ca, and 
portions with negative differential conductivity. As fol- 
lows from Eqs. (|4j) and (|55|), there are either one or two 
resonance voltages, depending on the relation between H 
and ui. The increase of the magnetic field broadens the 
peaks in J(V) which eventually disappear at higher H, 
because of the increase of the effective damping constant 
T] in Eqs. 

B. Ac dissipation 

The mean ac power Q = (JV) dissipated per unit area 
of a GB due to the ac voltage V(t) = V a cosuot of large 
amplitude can be calculated by solving Eqs. (|3l])-(|3^) 
numerically. The situation simplifies for a low-frequency 
ac signal (ujt <C 1) for which Q can be obtained using 
the quasistatic Eq. (|3^). For low fields ft<Cl and mod- 
erate ac amplitudes u a = V a /Vo < 1, the unity under the 
square root in Eq. d36j) can be neglected, giving 



Q = 



2irR 



2-n/u 



cos 2 Lotdt 



\Jh + u 2 cos 2 uot 



(65) 
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This integral can be expressed inj-terms of the complete 
elliptic integrals K(m) and E(m)cB: 



Q = 



:[(l+g 2 )E(m)-K(m)}, 



9 2 /(l+9 2 ), 



V a /R f J b . 



(66) 
(67) 



The power Q(g) is a function of only one dimensionless 
parameter g which includes both the ac amplitude and 
the magnetic field, as shown in Fig. 7. For a weak ac 
signal, V a <C JbRy/h, Eq. ( |66| ) yields the quadratic de- 
pendence of Q = V a 2 /2i?/ on V a , where Rf = R\fh is the 
dc flux flow resistivity at h <C 1. However, for stronger 
ac signals, JbRy/h < V„< J D R, Eqs. ( |65|) and ( |66| ) yield 
the linear dependence 



2V a J b /7T, 



(68) 



which is independent of the GB resistance R. This be- 
havior is due to the AJ core expansion at large vortex 
velocities. For very high ac amplitude V a 3> RJb, the 
full Eq. ( [36| ) should be used instead of Eq. j6q). In 
this case Eq. ( |36| ) yields V — J/R during most part of 
the ac cycle, thus the ac power dissipated on the GB 
becomes quadratic in V a , approaching the normal state 
limit, Q -> V a 2 /2R (not shown in Fig. 7). 



C. Generation of higher harmonics 

The nonlinearity of the electromagnetic response of 
the A J structure driven by harmonic ac current J(t) = 
J a cos Lot gives rise to higher voltage harmonics 



V{t) = ^V 2n +i{P,u,h) cos(2n+ l)wi, 



(69) 



where the Fourier coefficients V2 n +i can in principle be 
calculated from Eqs. (^0|)-(^2|). The higher harmonics 
in V(t) are most pronounced if the amplitude of ac AJ 
vortex displacements is maximum. For the overdampcd 
dynamics considered in this paper, the amplitude of the 
A J vortex oscillations decreases with u>, so the higher 
harmonics in V(t) are most pronounced for the quasi- 
static ac signal, lot C 1. In that case the coefficients, 

2lo 

V 2 n+i = — / dtV(/3 a cos Lot) cos(2n+ l)wt (70) 
it Jo 



are independent of lo, and V[/3(t)] is given by Eq. (35). 
For instance, the amplitude of the third harmonics V3 
for a small ac current f3 a — J a /Jb *C 1 is obtained by 
expanding Eq. ( |35| ) up to cubic terms in (3. This yields 
V(t) w RfJ a cosLut + V 3 cos3Lot, where 



V, = 



H 2 (i 



(71) 



The amplitude of the third harmonics V${h) has a max- 
imum at h = 1/4, as shown in Fig. 8. This field de- 
pendence reflects the increase of V3 oc y/h at small h 
proportional to the AJ vortex density, followed by the 
decrease of Vs(h) at higher fields, for which the AJ cores 
start overlapping, and the V(J) curve becomes Ohmic. 
This trend is characteristic of other higher order harmon- 
ics as well, whose amplitudes V2 n +i(h) oc f3 2n+1 strongly 
decrease with n and h for h > 1. 

The maximum in V^(h) could be used for extracting 
the field Hq and thus the AJ core size I on a GB from 
the ac measurements. In that case the ac measurements 
of V${H) may bring some advantages over the dc mea- 
surements of Hq from the flux flow resistivity Rf(H)c3, 
since the maximum in V^{H) occurs at the field Hq/A 
independent of the quasi-particle resistivity R, while the 
extracting of / from Rf(H) requires a two-parameter fit 
for Hq and R. 



D. Phase locking and quasi-steps on J — V curves 

AJ vortices in superimposed ac and dc fields can ex- 
hibit phase locking effects at ac large amplitudes (3 a if the 
dc voltage V on a GB is commensurate to the Joseph- 
son voltage V u = hio/2e, where n is any integer, and 
e is the electron charge. The electromagnetic response 
of a GB biased by a dc voltage V superimposed on ac 
voltage V a cosLot, is described by Eqs. (^l|)-(|33|) with 
u{t) = u+u a cosLot, where u = V/RJb, and u a = V a /RJb- 
To calculate the dc J — V characteristic averaged over 
the ac oscillations with the account of the phase locking 
of "/(t) onto the ac field, we use an approach similar to 
that for Shapiro cuccept steps on J — V curves of small 
Josephson junctionsHE3. In this case 

7 = ip + nuJ t + qsiriLot + 5j, q — u a /fl, (72) 

where ip is a constant phase shift, and S"f(t) is a nongrow- 
ing oscillating correction. Likewise, a = ao+da(t), where 
ctQ is a constant to be determined, and 5a(t) is an oscil- 
lating correction. We consider here the high-frequency 
signals of large amplitude (lot ^> 1, u a > 1), for which 
both 8a ~ 1/lot <C 1 and £7 ~ 1/lot -C 1 can be ne- 
glected. This situation differs from the nonlocked state 
(n = 0) considered in the previous section, for which the 
contribution of Sa and £7 entirely determine the effect of 
the ac field on the dc V — J curves. 

Eqs. @ and Eqs. @-@ yield the following aver- 
aged dc equations, 



(cos 7) sinh a>o = Vh, 
(sin 7) sinhao = u — nO, 
f3o = u+(sm 1 )e- ao . 



(73) 
(74) 
(75) 



Here the averages (sin 7) and (cos 7) were calculated in 
Appendix E: 



(sin 7) = (-l) n J n (q) sin^, 
(00S7) = (-l)"J„(g)cosV>, 



(76) 
(77) 
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where J n (q) is the Bessel function. From Eqs. (|73|)-(77) 
it follows that 



tan %p = (u- nQ)/Vh, (78) 
sinh 2 a Q = [h + (u- nQ.) 2 } / 'J%{q) . (79) 



Eqs. (|75|), (^8|), and (|79j) give the following dc J - V 
characteristics at V « nV u : 



RJ = nV u + (V - nV u ) 



h+(V -nv^y/v 2 



1/2 

(80) 



The behavior of J(V) near the resonant voltage V « nV^ 
is shown in Fig. 9. For H — ► Ho, the J(V) dependence 
approaches that of the Shapiro step for a small Josephson 
junction at zero field. For finite H , the contribution from 
the A J vortex motion broadens the kink in J(V), whose 
width AV = V — Vu, can be estimated from the condition 
that two terms in the denominator of Eq. ( p0| ) become 
comparable. This yields 



AV 



8ttA 2 ' 



(81) 



The value AV is independent of Jb and goes to zero as 
T c - T at T c . 

The evolution of the A J core length L = laj\fh as 
the dc voltage V sweeps through the resonance at V u is 
described by Eq. (|79|) as follows 



L 



I . i \\/h + (u — nil) 

—= sinh — : — 

Vh L \J n (q)\ 



2 , 



(82) 



The core length L(V) passes through minima at the res- 
onance voltages nHu>/2e. Notice that there are specific 
values of the ac parameter q n = u a /Q which correspond 
to zeros of the Bessel function J n (q) . In this case the core 
length diverges logarithmically, so the above approxima- 
tion, which neglects the ac corrections ~ 8a oc l/f2, 
becomes invalid. The blowing up of the AJ core length 
above A at q ~ * q m may indicate a conversion of the AJ 
into the J vortex under the action of an ac field. 



VII. DISCUSSION 

In this paper solutions that describe dc and ac driven 
mixed Abrikosov vortices with Josephson cores on high- 
Jb grain boundaries in a magnetic field are obtained. 
These solutions give self-consistent distributions of cur- 
rents circulating around moving AJ vortex structure in 
an exactly solvable model of the overdamped AJ vortex 
dynamics that describes both nonlinear dissipative pro- 
cesses in the vortex cores and magnetic interaction be- 
tween AJ vortices along a GB. Unlike Josephson vortices 
whose overdamped nonlinear dynamics in the long junc- 
tions can be described only numerically, the dynamics of 
AJ vortices turned out to be integrable just in the over- 
damped limit which is most relevant to low-angle GBs. 



The analytic theory of the ac response developed in this 
paper could be used to describe high- Jb flux flow oscilla- 
tors based on HTS bicrystalsH. 

Based on the exact A J dynamic solutions, both the 
dc flux flow resistivity and the V-J characteristics are 
obtained. The field dependence of the flux flow resis- 
tivity Rf{H) shows the characteristic \H behavior at 
low H, but then it approaches the quasiparticle limit R 
for H S> Hq as the A J cores overlap. The simplicity 
of Eq. ([s?]) for Rf(H) gives a direct way of extract- 
ing the AJ core length I and thus the intrinsic depair- 
ing current density Jb and the quasiparticle resistance 
R averaged over few current channels from transport 
measurements. Such measurements have indeed proven 
the existence of AJ vortices on 7° irradiated and unir- 
radiated YBCO bicrystals for which the AJ core lengths 
~ 100 — 200A at 55 — 77Kja considerably greater than 
£(T), but smaller than A(T)E3. In addition, the extracted 
temperature dependence of Jd(T) exhibited a clear SNS 
behavior Jb — Jo(l — T/T c ) 2 , indicating a significant 
order parameter suppression between dislocation cores, 
even on a rather low- angle 7° GB, in accordance with 
the model of Ref. ||. 

The fact that moving AJ vortex core can effectively 
probe local GB properties at the nanoscale of few GB 
dislocation spacings, makes standard transport measure- 
ments a very useful tool to clarify the structure of vortex 
core on GBs, mechanisms of current transport through 
GBs and the effect of local overdoping on Jb and R. 
This method implies measurements of flux flow resis- 
tance Rf(H) for bicrystals with the same misorienta- 
tion angle but different local doping level. In this 
way, the intrinsic depairing current density of a GB, Jb 
can be obtained as a function of the dopant concentra- 
tion. Pinning of AJ vortices strongly affects the V — J 
curve near the depinning current density, J ~ J g b, but 
for J 3> Jgb, the differential resistivity dV/dJ approaches 
the free flux flow resistivitjC3. Thus, measuring Rf(H) 
in the flux flow region at J > (2 — 3)J g b enables one 
to avoid the analysis of multiply pinning mechanisms on 
real grain boundaries in HTSB, using instead an exact 
flux flow theory for the interpretation of the experimen- 
tal data in the region J 3> J g b, where pinning is a weak 
perturbation. This conclusion is consistent with the ex- 
perimental fact that that V — J curves observed on HTS 
bicrystals ar.e r.athfjr straight above J g b in a wide range 
of currents^Baol^J ESEZrH. Thus, the intrinsic proper- 
ties of grain boundaries can be extracted from the anal- 
ysis of the differential resistance (but not V-J curves) 
measured at J ^ J g b using our solution for Rf which 
neglects pinning. A similar approach was used to mea- 
sure the flux flow resistivity of pinned A_vortices driven 
by strong current pulses well above J c Lil. Since J g b is 
by 2-3 orders of magnitudes below the intrinsic Jb, the 
pinning region is much smaller than the scale of Fig. 2, 
so the linear Rf can be used to fit the data. 

Measuring the local ac response of a GB can bring ad- 
ditional advantages over the dc measurements in which 
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the dc current acts both on AJ on GB and A vortices in 
the grains. Indeed, using a scanning localized microwave 
sourcelia, it would be possible to probe the ac dynamics 
of AJ on GB, minimizing the effect of pinned A vortices. 
For instance, the AJ core size and thus the local Jb can 
be extracted from ac measurements of higher voltage har- 
monics, as described in the previous section. 

The behavior of AJ vortices in a strong dc field H ^> 
H c i considered in this paper is most relevant for bulk 
HTS. However, of much interest for superconducting elec- 
tronics is also the electromagnetic response of HTS poly- 
crystaline films in comparatively weak rf fields for which 
vortices mostly penetrate the network of GBs. In equilib- 
rium, this regime corresponds to the field range H c ib < 
H <H cl , where H cX = (0 o /4ttA 2 )(1 - A)[ln(A/£W- 0.5] 
is the lower critical field for intragrain A vorticesEj, and 
H clb = (<?(> /47rA 2 )(l - N^n(X/l) + 0.423] is the lower 
critical field of AJ vorticeala, and N is the demagnetizing 
factor. For a film in a perpendicular field, the AJ vortex 
spacing a = (cpo/B) 1 / 2 can now significantly vary along 
GB, since the normal component of the local magnetic 
induction B(x,y) is now determined by highly nonuni- 
form distribution of the Meissner surface currents and 
a dome-like vortex distribution due to the geometrical 
barrierOEi In this case the results of this paper based 
on the solution (||) with constant k — 2n/a(B) may be 
used if B(x) varies slowly over the inter vortex spacing 
a(x), thus the parameter h in the above formulas should 
be replaced by its local value h(x) = B(x)/Hq. 
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APPENDIX A: DERIVATION OF EQ. (@) 

To calculate the 2D current distribution around a pla- 
nar Josephson contact in the xz plane at y = 0, it is 
convenient to use the scalar stream function ip(x,y) so 
that J x = d y ip, J y — —d x ip, whence 



4ttA 2 



dip 
dy 

dtp c 
~dx = ~4vrA 2 



r I 0Q dp 

2tt dx 

(p dp 



2tt dy Ay 



(Al) 
(A2) 



where A is the vector potenti al, and p is the phase of the 
order parameter. From Eq. (Al), it follows that the par- 
allel component J x (x,+0) — J x (x,—0) = (c<p / Sir 2 \ 2 )6' 
is discontinuous for any nonuniform distribution of the 
phase difference 9{x) — ip(x, +0) — p(x, —0) along GB. 
This results in the following boundary condition for the 
stream function at a GB: 



d y i/>(x, +0) - d y iP(x, -0) = (c0o/87r 2 A 2 )6>' 



(A3) 



Excluding p from Eqs. ( |Al[ ) and (A2), we obtain 
cH 



4ttA 2 ~ 8tt 2 A 2 



(A4) 



where the delta function 8(y) in the right-hand side pro- 
vides the boundary condition (A3), and H = V z x A 
is the z-component of the magnetic field, related to the 
stream function ip by the Biot-Savart law: 



H(r) = \ 



{(x-x>)d x ,yj + (y~y')d y ^}d 3 r' 
[(x - x') 2 + (y- y') 2 + (z- z') 2 } 3 / 2 ' 



(A5) 



Eqs. (A4) and (A5) give an integro-differential equation 
for ip(x, y) which can be solved by the Fourier transform. 

For an infinite (along H) sample, the local relation 
ip(x,y) = cH(x,y)/4:TT holds, thus Eq. (A4) becomes 
the London equation ( |l5| ) for H(x,y). For a thin film 
of thickness d <C A, one can put z = z' = in Eq. 
(|AE| ) , then the integration over z' gives the factor d, and 
Eq. (A5) gives the following relation between the Fourier 
components of H and ip: 



(2tt qd/c)ip q , 



(A6) 



where q 2 
yields 



q . The Fourier transform of Eq. (A4) 



ccp Q 6'(q x 



8ir 2 \ 2 (q 2 +q/X) 



(A7) 



where A = 2 A 2 jd. Making the inverse Fourier transform 
of Eq. ([AT]) and using the boundary condition ip'{x, 0) = 
Jftsinf? + HG/2eR, we arrive at the nonlocal Eq. (|4|) for 
6(x, t) in which Kq(\x — wl^J^hould be replaced by the 
kernel K(\x — u\/X), where 



K(s) = 



-Sp 



dp 



(A8) 



Despite different behaviors of KqQx— u\/X) and K(\x— 
u\/X), at large distances, they both have the same loga- 
rithmic singularity at x = u, because at short distances 
\u — x\ <C A, scr eenin g is inessential so the vector poten- 
tial A in Eqs. jAl|)-(|A^) and the field H in Eq. @) 
can be neglected. Then Eq. (A4) becomes a 2D Poisson 
equation whose solution is 



ip{x,y) 



32tt 3 A 2 



9'(u)ln[y 2 + (x-u) 2 ]du (A9) 



both for thin films and bulk samples. 



APPENDIX B: EXACT SOLUTION 

To show that Eq. (^|) is indeed an exact solution of Eq. 
]), we substitute Eq. @-@ into Eq. (0) in which the 
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integral is evaluated using the Hilbert transfer: 
sinh ady ir sin kx 



_ QO (x — y) (cosh a — cosky) cosh a — cosfcx 



(Bl) 



Eq. (0) then reduces to the algebraic form C\ cos ( + 
C2 sin C + C3 = 0, where the coefficients Cj depend only 
on time via a(t), "f{tV and Xo(t). Equating all Cj to zero, 
we arrive at Eqs. ©-(jl|). 

Useful relations for the steady-state vortex propaga- 
tion can be obtained from Eqs. (10)-(12|) with 7 = a = 0: 



sinh a cos 7 = 
sin 7 cosh a = /3 (* ) , 
svh = — sin 7 sinh a. 



Adding squared Eqs. (|B2| ) and (B4) yields 

sinh 2 a = (1 + s 2 )h. 
Subtracting squared Eqs. ( p33| ) and (B4) yields 
sin 2 7 = 0q- s 2 h, 



(B2) 
(B3) 
(B4) 



(B5) 



(B6) 



Substituting Eqs. (|BJ) and back to Eq. (|BJ) gives 
the following equation for s 



/9g = /is 2 + 



1 



(B7) 



The use of the relation syh = V/ J^R reduces Eq. (B7) 
to Eq. d|). Substituting Eq. @ into Eq. ® yields 



sin 7 



cos 7 = 



vT 



(B8) 



These results can also be obtained from the complex 
representation of the dc equations p4| ) 

/9b + iVh. (B9) 

The vortex velocity is given by the second Eq. (jl4|) 

sVh = Im[l - (A, + iV77)] 1/2 , (BIO) 

from which Eq. fl34| ) readily follows. 

APPENDIX C: R f FOR J VORTICES 

The sine-Gordon equation for a periodic J vortex struc- 
ture moving with a constant velocity v has the form 



X 2 d C( 9 + vrd c 9 - sin 6> + (3 = 0, 



(CI) 



where the £ — x — vt. Multiplying Eq. (CI) by 8^9 
and then integrating from ( = to ( = a, using the 
flux quantization condition 9(C, + a) — 9(() = 2tt, and the 



periodicity condition d^9(( + a) = d^9(() for -ff (£), we 
arrive at the following equation for v: 



(C2) 



2tt/3 = -vt / (d c eydc 



Expressing the vortex velocity v(J) in Eq. (|C2j) via the 
dc voltage V from Eq. ( |30| ) yields the flux flow resistivity 
Rf = V/ J in the form 

R f = ATt 2 R/\a f (d ( 9) 2 d(\ (C3) 

For the linear flux flow state at f3 < 1, 6>(C) in Eq. (|c|) 
can be replaced with static 9q(x) for which 

AX'- sin O = O, ( C4 ) 
The first integral of Eq. ( |C4| ) has the form 

(Aj0 o ) 2 = 4[p- 2 - cos 2 9 /2] (C5) 



The periodic solutions of Eq. ( |C5| ) can be expressed 
in terms of the Jacobi elliptic function, 9q(x) = n + 
2am(x /p\j), where < p < 1 is a parameter which 
is related to the applied magnetic field H and the mag- 
netic induction B by Eqs. (|39|)E3. The period a of the J 
structure is given by 

a = 2\ jP K(p), 2\aB = tfo, (C6) 

where Kip) and E(p) are the complete elliptic 



integrals 



r /2 

K(p)= / (l-p 2 sm 2 9)- 1/2 d9, (C7) 
Jo 

E{p)= (1 -p 2 sin 6 2 ) l ' 2 d6 (C8) 
Jo 



Substituting Eq. ( |C5| ) into Eq. (C3) and changing inte- 
gration J dx — > J Q d9/9' , we obtain 

ir 2 Rp\j 



Rf = 



2aE{p) 



(C9) 



Eqs. flCfl), (C9) and (B9J) give Eq. (B 



APPENDIX D: LINEAR RESPONSE 



Linearized Eqs. (J10|)-(|l2|) have the form 

r<5ci + 8a cosh a cos 7 — ^7 sinh a sin 7 = 0, (Dl) 

T(57 + <57COS7Cosha + 5a; sin 7 sinh a = /3 w e lwt , (D2) 

krSv = —^7 cos 7 sinh a — 8a sin 7 cosh a (D3) 

These equations give the following Fourier components 
of 8a, £7, and <5s: 

sinh a sin 7<5/3 w 



2 a - ' ( D4 ) 

(if2 + cosh a cos 7) 2 + sinh a sin 7 



(id + cosh a cos 7)5/3^ 



VhSi 



(XI + cosh a cos 7) 2 + sinh 2 a sin 2 7 

sinh a(i£7 cos 7 4- cosh a)5/3 w 
(ifl + cosh a cos r j) 2 + sinh 2 a sin 2 7 



, (D5) 



, (D6) 
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where Q — ut. Substituting Ss and £7 into Eq. j30| ) 
gives the complex ac resistance — V u / J a : 



R 



sinh a cosh a — ft 2 + ifle a cos 7 
sinh 2 a + cos 2 7 — £1 2 + 2iQ cosh a cos 7 



(D7) 



Eq. (D7) reduces to Eq. (M), using Eqs. (B5)-(B8) 



APPENDIX E: SEPARATION OF FAST AND 
SLOW VARIABLES 

The time averages of 

(sin 7) = ( (sin ip cos mot + cos ip sin nut) cos(q sin tot) 

+(cos"0 cosnwi — sin'0 sinnuit) sin(g sinwt))(El) 

(cos 7) = ((cos ip cosnut — sin?/' smmot) cos(q shout) 

— (sin-0 cos mot + cos ip sinmot) sin(q sinwi))(E2) 



can be calculated, using the identities^ 



cos(q sin tot) = J (q) + 2 J 2 fe (9) cos(2fcwi) , (E3) 

00 

sin(qsmujt) = 2 J2fc+i(g) sin[(2fc + l)wt], (E4) 



fc=0 



where Jk{q) is the Bessel function. Averaging Eqs. (El) 
and (E2) yields Eqs. ( |76| ) and ( f77j ) because only one 
resonant term with 2fc = n or 2k + 1 = n in the sums 



and (cos 7). 



( E3 ) and ( E4 ) gives a nonzero dc contribution to (sin 7) 



1 C.C. Tsuei and J.R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

2 D. Larbalestier, A. Gurevich, D.M. Feldman, and A. 
Polyanskii, Nature 414, 368 (2001). 

3 P. Chaudhari, J. Mannhart, D. Dimos, C.C. Tsuei, J. Chi, 
M. Oprysko, and M. Scheuermann, Phys. Rev. Lett. 60, 
1653 (1988); D. Dimos, P. Chaudhari, and J. Mannhart, 
Phys. Rev. B 41, 4038 (1990). 

4 N.D. Browning, J. P. Buban, C. Plauteau, G. Duscher, 
and S.J. Pennycook, Micron 30, 425 (1999); S.E. Bab- 
cock, ibid., p. 449; A. Goyal, S.X. Ren, E.D. Specht, D.M. 
Kroeger, R. Feenstra, D. Norton, M. Paranthaman, D.F. 
Lee, and D.K. Christen, ibid., p. 463. 

5 K.L. Merkle, Y. Gao, and B.V. Vuchic, in Characteri- 
zation of High T c Materials & Devices by Electron Mi- 
croscopy, Edited by N.D. Browning and S. Pennycook 
(Cambridge University Press, Cambridge, New York, Mel- 
bourne, Madrid, 2000), p. 235; N.D. Browning, M.F. 
Chisholm, and S.J. Pennycook, ibid., p. 263; 

6 A.P. Sutton and R.W. Ballufh Interfaces in Crystalline 
Materials, (Clarendon, Oxford, 1995). 

7 H. Hilgenkamp and J. Mannhart, Rev. Mod. Phys. 74, 485 
(2002). 

8 S.E. Russek, D.K. Lathrop, and B. Moeckly, Appl. Phys. 
Lett. 57, 1155 (1990). 

9 A. Schmehl, B. Goetz, R.R. Schultz, C.W. Schneider, 
H. Bielefeldt, H. Hilgenkamp, and J. Mannhart, Euro- 
phys. Lett. 47, 110 (1999); G. Hammerl, A. Schmehl, 
R.R. Schultz, B. Goetz, H. Bielefeldt, C.W. Schneider, H. 
Hilgenkamp, and J. Mannhart, Nature 407, 162 (2000). 

10 G.A. Daniels, A. Gurevich, and D.C. Larbalestier, Appl. 
Phys. Lett. 77, 3251 (2000). 

11 K. Guth, H.U. Krebs, H.C. Freyhardt, and Ch. Jooss, 
Phys. Rev. B 64, R140508 (2001). 

12 A. Gurevich and E.A. Pashitskii, Phys. Rev. B 57, 13878 
(1998). 

13 J. Mannhart and H. Hilgenkamp, Mater. Sci. Engr. B 56, 
77 (1998); Physica C 317-318, 383 (1998). 



14 J. Halbritter, Phys. Rev. B 48, 9735 (1993); J. Supercond. 
10, 91 (1997); Phys. Rev. B 61, 1596 (2000). 

15 S.M. Anlage, W. Hu, CP. Vlahacos, D. Steinhauer, B.J. 
Feenstra, S.K. Dutta, A. Thanawalla, and F.C. Wellstood, 
J. Supercond. 12, 352 (1999). 

16 Y.M. Habib, C.J. Lehner, D.E. Oates, L.R. Vale, R.H. 
Ono, G. Dreesselhaus, and M.S. Dresselhaus, Phys. Rev. 
B 57, 13833 (1998). 

17 M. Golosovsky, M. Tsindlekht, and D. Davidov, Super- 
cond. Sci. Technol. 9, 1 (1996); M. I. Tsindlekht, E.B. 
Sonin, M.A. Golosovsky, D. Davidov, X. Castel, M. 
Guilloux-Viry, and A. Perrin, Phys. Rev. B 61, 1596 
(2000). 

18 V.P. Koschelets and S.V. Shitov, Supercond. Sci. Technol. 
13, R53 (2000). 

19 A. Gurevich, Phys. Rev. B 46, R3187 (1992); 48, 12857 

(1993) . 

20 A. Gurevich, Physica C243, 191 (1995). 

21 A. Gurevich and L.D. Cooley, Phys. Rev. B 50, 13363 

(1994) . 

22 A. Diaz, L. Mechin, P. Berghuis, and J.E. Evetts, Phys. 
Rev. Lett. 80, 3855 (1998); Phys. Rev. B 58, R2960 (1998). 

23 D.T. Verebelyi, D.K. Christen, R. Feenstra, C. Cantoni, 
A. Goyal, D.F. Lee, M. Paranthaman, P.N. Arendt, R.F. 
Depaula, J.R. Groves, and C. Prouteau, Appl. Phys. Lett. 
76, 1755 (2000); D.T. Verebelyi, C. Cantoni, J.D. Budai, 
D.K. Christen, H.J. Kim, and J.R. Thompson, ibid. 78, 
2031 (2001). 

24 D. Kim, P. Berghuis, M.B. Field, D.J. Miller, K.E. Gray, 
R. Feenstra, and D.K. Christen, Phys. Rev. B 62, 12505 
(2000). 

25 J. Albrecht, S. Leonardt, and H. Kronmuller, Phys. Rev. 
B 63, 014507 (2001). 

26 M.J. Hogg, F. Kahlmann, E.J. Tarte, Z.H. Barber, and 
J.E. Evetts, Appl. Phys. Lett. 78, 1433 (2001). 

27 H. Claus, U. Welp, H. Zheng, L. Chen, A.P. Paulikas, B.W. 
Veal, K.E. Gray, and G.W. Crabtree, Phys. Rev. B 64, 
144507 (2001). ' 



15 



28 E.L. Kramer, J. Appl. Phys. 41, 621 (1970). 

29 H.J. Jensen, Y. Brychet, and A. Brass, J. Low Temp. Phys. 
74, 293 (1989); H.J. Jensen, A. Brass, A.-C. Shi, and A.J. 
Berlinski, Phys. Rev. B 41, 6394 (1990). 

30 C. Reichhardt, C.J. Olson, and F. Nori, Phys. Rev. B 57, 
7937 (1998); 58, 6534 (1998); C. Reichhardt, G.T. Zi- 
manyi, and N. Gronbech- Jensen, Phys. Rev. B 64, 015401 
(2001). 

31 G.W. Crabtree, D.O. Gunter, H.G. Kaper, A.E. Koshelev, 

G. K. Leaf, and V.M. Vinokur, Phys. Rev. B 61, 1446 
(2000). 

32 C.A. Bolle, PL. Gammel, D.G. Grier, C.A. Murray, D.J. 
Bishop, D.B. Mitzi, and A. Kapitulnik, Phys. Rev. Lett. 
66, 112 (1991). 

33 T. Matsuda, O. Kamimura, H. Kasai, K. Harada, T. 
Yoshida, T. Akashi, A. Tonomura, Y. Nakayama, J. Shi- 
moyama, K. Kishio, T. Hanaguri, and K. Kitazawa, Sci- 
ence 294, 2136 (2001). 

34 A. Pruijmboom, P.H. Kes, E. van der Drift, and S. Rade- 
laar, Phys. Rev. Lett. 60, 1430 (1988). 

35 R. Besseling, R. Niggerbrugge, and P.H. Kes, Phys. Rev. 
Lett. 82 3144 (1999); S. Anders, A.W. Smith, R. Besseling, 
P.H. Kes, and H.M. Jaeger, Phys. Rev. B 62, 15195 (2000). 

36 A. Gurevich, Phys. Rev. B 42, R4857 (1990); A. Gurevich, 

H. Kiipfer and C. Keller, Europhys. Lett. 15, 789 (1991). 

37 R. Wordenweber, Phys. Rev. B 46, 3076 (1992). 

38 M. Ziese, Phys. Rev. B 53, 12422 (1996); 55, 8106 (1997). 

39 K. Yamafuji and T. Kiss, Physica C258, 197 (1996). 

40 M. Prester, Supercond. Sci. Technol. 11, 333 (1998). 

41 M. Suenaga, in Superconductor Materials Science, edited 
by S. Foner and B.B. Schwartz (Plenum, New York, 1981), 
p. 201. 

42 E.L. Kramer and H.C. Freyhardt, J. Appl. Phys. 51, 4930 
(1980); C. Rossel and O. Fisher, J. Phys. F14, 455 (1984). 

43 A. Barone and G. Paterno, Physics and Applications of the 
Josephson Effect (Wiley, New York, 1982). 

44 K. K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

45 Yu.M. Ivanchenko and T.K. Soboleva, Phys. Lett. A147, 
65 (1990). 

46 Yu.M. Aliev, V.P.Silin, and S.A. Uryupin, Superconduc- 
tivity, 5, 230 (1992); JETP Lett. 57, 193, (1993); Yu.M. 
Aliev and V.P. Silin, JETP 77, 142, (1993). 

47 R.G. Mints and LB. Snapiro, Phys. Rev. B 51, 3054 (1995). 

48 V.G. Kogan, V.V. Dobrovitski, J.R. Clem, Y. Mawatari, 
and R.G. Mints, Phys. Rev. B 63, 144501 (2001). 

49 G. L. Alfimov and A. F. Popkov, Phys. Rev. B 52, 4503 
(1995). 

50 Yu.E. Kuzovlev and A.I. Lomtev, JETP 84, 986 (1997). 

51 A. Seeger, Theorie der Gitterfehlstellen, In Handbuch der 
Physik (Springer, Berlin, 1955), vol. 7, part 1. 

52 G.L. Alfimov and V.P. Silin, JETP 79, 369 (1994). Note 
added: after this manuscript was submitted I became 
aware of the papers by V.P. Silin, JETP Lett. 60, 462 
(1994) and JETP 85^ 760 (1997) in which equations anal- 
ogous to Eqs. @- ( |l2] ) and ( |35[ ) had also been obtained. 

53 A. Gurevich, M.S. Rzchowski, G. Daniels, S. Patnaik, B.M. 
Hinaus, F. Garillo, F. Tafuri, and D.C. Larbalastier, Phys. 
Rev. Lett. (2002). 

54 J. Bardeen and M. J. Stephen, Phys. Rev. 140, 1197 
(1965). 

55 A.T. Fiore, Phys. Rev. Lett. 27, 501 (1971); Phys. Rev. B 
7, 18811 (1973); 8, 5039 (1973). 

56 P. Martinoli, O. Daldini, C. Leemann, and B. Van den 
Brandt, Phys. Rev. Lett. 36, 382 (1976); P. Martinoli, 



Phys. Rev. B 17, 1175 (1978). 

57 R.E. Eck, D.J. Scalapino, and B.W. Taylor, Phys. Rev. 
Lett. 13, 15 (1964). 

58 M. Girillo, M. Gronbech- Jensen, M.R. Samuelsen, M. 
Salerno, and G. Verona Rinati, Phys. Rev. B 58, 12377 

(1998) ; M. Salerno and M.R. Samuelsen, ibid. 59, 14653 

(1999) . 

59 G.E. Blonder, M. Tinkham, and T.M. Klapwijk, Phys. 
Rev. B 7, 4515 (1982). 

60 V.L. Pokrovsky and A.L. Talapov, JETP 48, 971 (1978); 
51, 134 (1980). 

61 LM. Floria and J.J. Mazo, Adv. Phys. 45, 505 (1996). 

62 X.Y. Cai, A. Gurevich, I-Fei Tsu, D.L. Kaiser, S.E. Bab- 
cock, and D.C. Larbalestier, Phys. Rev. B 57, 10951 
(1998). 

63 F.F. Ternovskii and L.L. Shekhata, Zh. Exp. Teor. Fiz. 
62, 2292 (1972) [Sov. Phys. JETP 35, 1202 (1972)]; A.E. 
Koshelev, Physica C223, 2761 (1994). 

64 J.I. Gittleman and B. Rosenblum, Phys. Rev. Lett. 16, 734 
(1966). 

65 M.W. Coffey and J.R. Clem, Phys. Rev. Lett. 67, 386 
(1992); J.R. Clem and M.W. Coffey, Phys. Rev. B 46, 
14662 (1992). 

66 J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 

67 P. Lebwohl and M. J. Stephen, Phys. Rev. 163, 376 (1967); 
J.R. Waldram, A.B. Pippard, and J. Clarke, Phil. Trans. 
Roy. Soc. (London) A 268, 265 (1970). 

68 I.O. Kulik, Zh. Exp. Teor. Fiz. 51, 1952 (1966); C.S. Owen 
and D.J. Scalapino Phys. Rev. 164, 538 (1967). 

69 Handbook of Mathematical Functions., Applied Mathemat- 
ics Series, vol. 55, Edited by M. Abramovitz and LA. Ste- 
gun (National Bureau of Standards, 1964). 

70 G. Blatter, M.V. Feigelman, V.B. Geshkenbein, A.I. 
Larkin, and V.M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 

71 M.N. Kunchur, D.K. Christen, and J.M. Phillips, Phys. 
Rev. Lett. 70, 998 (1993); M.N. Kunchur, B.I. Ivlev, D.K. 
Christen, and J.M. Phillips, Phys. Rev. Lett. 84, 5204 

(2000) . 

72 E.H. Brandt and M.V. Indenbom, Phys. Rev. B 48, 12893 
(1994); 

73 E. Zeldov, J.R. Clem, M. McElfresh, and M. Darwin, Phys. 
Rev. B 49, 9802 (1994); E. Zeldov, A.I. Larkin, V.B. 
Geshkenbein, M. Konczykowski, D. Mayer, B. Khaykovich, 
V.M. Vinokur, and H. Shtrikman, Phys. Rev. Lett. 73, 
1428 (1994). 



16 



FIG. 1: Current streamlines around A J vortices on a GB 
(dashed line) and the bulk A vortices in the grains, calculated 
from Eq. (fjj) for I = 0.2a. 

FIG. 2: The V — J curves calculated from Eq. (||) for dif- 
ferent magnetic fields h = H/H : 0.01(1), 0.05(2), 0.1(3), 
0.5(4), 10(5). Inset shows the field dependence of the flux 
flow resistance Rf(B) given by Eq. ([jt]). 



FIG. 3: Dependencies of Qo and r\ on the normalized dc cur- 
rent density /3q = J/Jb for H = 0.05i?o- 
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FIG. 4: Field dependence of the ac dissipated power Q(H) 
for different dimensionless frequencies Q — lot. 



FIG. 5: Flux flow resonance line in Imii w (fi) calculated from 
Eq. @ for ft = 4 and different magnetic fields H/H : 0.05, 
0.4, and 2 (from bottom to top curve, respectively). 



FIG. 6: Manifestation of the flux flow resonance on the aver- 
aged dc J — V characteristic for cut — 0.1, H — 0.01/fo, and 
different ac amplitudes, (J a /2J b ) 2 : (1); 0.05 (2); 0.1 (3). 



FIG. 7: The power Q dissipated on a GB as a function of the 
ac amplitude, where Q a = 2J^Rf/n. 



FIG. 8: The field dependencies of the amplitude of the third 
harmonics Vs(H) normalized to V n = RJt/l&Jl and the flux 
flow resistance Rf(H). 



FIG. 9: Quasi-step on the averaged J ~ V characteristic at 
V fa Vu> for q = 3, n = 1, Vu,/RJb = 1, and different magnetic 
fields H/H : (1), 0.01 (2), and 0.1 (3). 
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